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摘 E 本 文 研究 了 二 维 空间 中 光滑 的 凸 区 域 上 具有 小 周期 参数 的 抛物 型 方程 的 半 离 散 双 义 度 有 限 元 坎 
近 。 在 整个 区 域 上 ， 利 用 二 次 元 解 均匀 化 问题 ， 在 参考 单 胞 内 ， 利 用 线性 元 解 辅助 问题 ， 然 后 ， 
将 求 得 的 有 限 元 解 代 入 多 尺度 渐 近 展开 式 ， 得 到 原 问题 解 的 一 个 半 离 散 双 尺度 有 限 元 格式 。 利 用 
多 尺度 渐 近 展开 和 有 限 元 理论 ， 证 明了 该 格式 的 收敛 性 。 
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1 引言 


随 着 计算 机 软 硬 件 的 发 展 ， 出 现 了 许多 涉及 具体 物理 问题 的 算法 与 理论 。 但 对 于 小 周期 复 
合 材 料 的 计算 ， 车 0 < = « 1 时 ，a$;(z) 变化 非常 频繁 ， 用 差分 或 有 限 元 等 数值 求解 时 ， 网 格 
前 分 都 要 求 非常 细 ， 这 样 导 臻 最终 求 解 线性 方程 组 的 计算 量 非常 大 。 针对 此 类 问题 ， 发 展 起 来 
了 均匀 化 和 多 尺度 渐 近 展开 方法 ， 即 在 宏观 尺度 上 求 均匀 化 问题 ， 在 微观 尺度 上 增加 校正 项 ， 
用 均匀 化 解 和 校正 项 构成 原始 解 的 渐 近 展开 式 ， 从 而 降低 了 计算 量 。 OCEA [1] 中 给 出 了 均匀 化 
理论 并 对 椭圆 型 方程 作 了 多 尺度 渐 近 展开 式 ， 文 献 [2-7] 从 不 同方 面 对 文献 [1] 中 的 结果 作 了 改 
进 。 此 外 ， 文 献 [8] 中 又 给 出 了 解决 椭圆 型 方程 的 一 种 新 思路 ， 文 献 [9,10] 又 利用 HMM 方法 分 
别 讨 论 了 抛物 型 方程 和 Helmholtz 方 程 的 多 尺度 计算 。 文 献 [11] 中 首次 给 出 抛物 型 方程 的 渐 近 
展开 式 和 收敛 性 分 析 。 本 文 利用 文献 [2] 中 提出 的 双 尺 度 有 限 元 方法 和 文献 [8-13] 中 的 思想 ， 对 
抛物 型 方程 给 出 一 种 半 离 散 双 尺度 有 限 元 格式 及 其 收敛 性 分 析 。 


2 ”预备 知识 和 已 有 结论 


本 文 考虑 如 下 问题 
u(x,t) — V - (at(z) Vue(z,t)) = f(z,t), (z,t) eO x (0,7), 
u, (2, t) = 0, (x,t) € 0€ x (0, T), (1) 
ug(z,0) = uł (zx), TEQ. 
其 中 
usto MED sr) aila 49 = 1,2, 
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它们 是 以 < 为 周期 的 周期 函数 。 设 y = 三， 车 (aug(y))2x2 € M(a, 8) 表示 存在 常数 0 < a < 
p. TER 
2 2 
oM n? Sag(y)yun; S 8 PR, V (m.m) # (0.0). (2) 
il i=1 


文中 采用 Einstein 求 和 记号 ， 出 现 的 Q = [0,1] x [0,1] 均 表示 参考 单 胞 ， 常 数 C 大 小 可 能 不 
相等 ， 但 都 与 = 无 关 。 为 使 考虑 的 问题 是 适 定 的 ， 设 六 L'(0,7; H^ (Q)). u(x) e LN) B. 
"ix € ON u(x) = 0。 此 外 ，(.,)p ERED) PINAR, |l- p KERHD) 中 的 模 ， 当 下 
标 为 只 时 可 以 省 略 。 

由 均匀 化 理论 知 ，(1) 的 均匀 化 问题 为 

u (x,t) - V- (a9 Vu(z,t)) = f(z,t), (x,t) € 0 x (0,7), 
u(x,t) = 0, (x,t) € 90 x (0, T), (3) 


u(a,0) = u(x), ce, 


这 里 pn 

, 
m i a0 = (a2;)2x2; aj; = E (a; — Qik ay v 

HP Ne, Nue HEQ) 满足 

a ie Oaik 





Oy; ne Oy; Oy 
0 ONxu 0 O(aij9k; Ni) m 
~ Oy; (a 5 Oy; =) = -ap - ae a Oy; Ye (4) 
Buls, t) AA FRÜUEGNTBDEOROTON 
Ou(z, u(x,t 
ul (a,t) = u(z,t) — ENk et +e ee. (5) 


利用 文献 [11] 的 证 明 方法 ， 可 得 : 
引 理 2.1 设 Q 是 R? 中 的 光滑 凸 区 域 


feCc(0T;H*(Q), as(z) EL™(0), üj-12, 


H (a£; (z))2x2 € M(a, 8). ue(2,t) 渐 近 展开 式 如 (5) 所 示 。 如 果 e CH(0, T; H4(Q)), Ne, New 
€ Wilo(Q) (kd = 2)， 则 有 下 面 的 估计 


< Ce CEU, (6) 


L2(0,T;H1(0)) — 





Pule, DM an, Pulz, t) 
DE Na DZKOT; ) 








uc(z,t) — (ule, t) - EN. 








这 里 C 与 = 无 关 。 
利用 文献 人] 中 引 理 3.2 的 证 明 方法 ， 可 以 证 明 : 
引 理 2.2 ||N.(Z)]lon < Ciel Nx (£)h;o < C2 B.1Nk(2)l20 S Ce-2, k 21,2. 
引 理 2.3 由 Nuez(E)llo， as Cie| Ni( 2 Vine € Co H |N (E )lan € CET 2? k,l— 13,2. 
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3 ”高 - 低 阶 的 双 尺 度 有 限 元 分 析 


由 渐 近 展开 式 (5) 知 ， 若 要 获得 问题 (1) WARR RARA RERE ulr, t), Ne 和 Nx。 
KERMA Lagrange WEN ERKE ulz, t), MAARTE Lagrange TE Q ERAR Ny 和 Nu 。 

KT QI Sc. MIERA A BV. CHQ ET, EYED RE RGE. AE 
4 Ny 的 有 限 元 近似 为 : SK NP € Vs, 满足 


(mE. Sea = (aix, Seo Vue Vp. (7) 


利用 经 典 的 有 限 元 证 明 方法 ， 可 以 得 到 下 面 结果 : 
引 理 3.1 ÆN, e W2(Q)， 则 有 


INe — Nellz2(Q) + hI Ne — Nea) < Ch?INi Emo); 
IN — Ne llr=(Q) + ALN ~ NE |wa(qy € Ch? In hINk|wz Q 


为 近似 求解 均匀 化 问题 (3)， 我 们 用 NT 代替 
0 _ 0 Ü s T3 ON; 
a = (aij)2x2 a; = Jo aa rh )dy 


中 的 Nj 得 到 a? 的 一 个 近似 a8 = (a5), itr 


ON? 
O,h 了 
Ql 一 一 Qij — Qik dy. 
1j ji ( J T Yk ) 


利用 文献 [2] 中 的 方法 ， 可 以 证 明 下 面 结论 : 
引 理 3.2 KN; e H?(Q), WA 


la; — ad| € Ch? Ni qo) N| is (g- 
VO RS — S 807A TA ho ARDER, BV, Ah E R3 — CUR DE I 26 TRI .ONL 
于 (3) 式 ， 其 半 离 散 双 尺度 有 限 元 格式 为 


(Bh (2, t), v) + (a2 Vün,, Vv) = (fov), Vv € Vig: 
W (8) 
T(x, 0) = uf, (x), V rcf. 
u(x) € H3(Q) B. 
Vu? (£) ~ Vuh, (z)lloa € Challu°(z)II3,0, 
lu? (x) — uho (z)]lo. < Chl|u?(z)]a.o; 


利用 引 理 3.1， 引 理 3.2 和 文献 [12] 中 的 证 明 技巧 ， 则 有 下 面 结果 : 
引 理 3.3  Wu(z, t) M an (x, t) PHA (3) M (8) E, F ulz, t) € H3(Q), Ns e H?(Q). 
则 有 


lu Ss. Un, || rz (o, (9) < C(h? + h2). 
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ON Ov — on —O(aió NP) O(NP — y) 
(u A E ( any 一 By. aki 2m o) VocV,. (9) 


利用 标准 的 有 限 元 理论 ， 引 理 3.1 和 引 理 3.2， 则 有 下 面 结果 。 
引 理 3.4 KN, e H?(Q)(k = 12)， 则 有 


[Ne — Nille € ChINI;o; 
Nav — Néllwa o) € Chll Nillwa Q). 


最 后 ， 我 们 给 出 ve 的 近似 计算 公式 为 








Folre HT), TET,» EX 


lloll no = "a 5 hnd)’. 


TETho 


定理 3.1 设 u(z,t) € C' (0,7; WE (9) n H* (0) 28 RUBE (3) HO, (x,t) € H?(0) 23 (8) XX 
的 解 。 # h < ho, Nx(y); Nxi(y) € w2 (Q) (k,l = 1,2), 则 
Juelz, t) — ug lin, < Cle? +h + h2). 


证 明 我们 有 


lue Cz, t) — us lli 


li 


h 
jue (ar, t) — uz (x,t) + uz (a, t) — ug llo 























< |jue(z, t) — uk (a, DN,no + ud (a, t) — us It ho 

< Ce + |lu — Gy, llz2(o rmi (o) + sw. - NR en TNT 
«eve (zs - 5) bus ae "ave - Nba L2(0,T;H}(Q)) 
«ela Lx 2522] | Iho 


{> 


= Ce? + A + A2 + 3 + A4 As. 


对 上 面 定义 的 A, Ao, As, Aa, 45， 可 分 别 估计 如 下 : 
对 于 4:， 由 引 理 3.3 可 得 4 < C(h? + FG). 
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WF A. A 


Ag 


IA 


Ce max, (Ni — Nelualulzsq,rwa, c) + Ne — Ne lolulz2@,7,w2,(a) ) 


IA 


Ce max, (e! Nx — Niels olulzz(or wa, co + Ne — NElz2qylelz2@,7,w2,a) } 


lA 


C max, {hl Nel2,qlelz2e,r;w2,a)) + Ek? 1Nx T2, lil zo, w2,()) 


lA 


C(h - eh?) € C(h - & 4 h*) € C(h - e). 

同样 ， 对 于 As, Aa, As 分 别 有 

Aa € C(eho + h2 +h?) < C(e? +h? + h?), 

Ag € C(?h? + eh) € Cet + ht +e? +h?) < Cle? +h?), As < Ce”. 


Se EATS, Juelz t) — ug” lias € C(e3 +h + h2)- 证 毕 


0 = 
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Two-scale FEM Analysis for Parabolic Equations with 
Small Periodic Coefficients 
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Abstract: We study the parabolic equation with small periodic coefficients in a smooth convex do- 
main of two-dimensional space and derive a semi-discrete two-scale finite element approximation in the 
paper. In the whole domain, we solve the homogenization problem with a quadratic element. In the 
reference cell, we solve the auxiliary problem with a linear element. Using the obtained finite element 
solutions to replace the solutions in the multi-scale asymptotic expansion, we get a semi-discrete two- 
scale finite element scheme for the original problem. Using the multi-scale asymptotic expansion and 
the finite element theory, we prove that the new scheme is convergen. 

Keywords: multiscale asymptotic expansion; two-scale FEM; parabolic equations; composite materi- 


als 
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